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The paper presents the group theory of best localized and symmetry-adapted Wannier functions 
in a crystal of any given space group G or magnetic group M. Provided that the calculated band 
structure of the considered material is given and that the symmetry of the Bloch functions at 
all the points of symmetry in the Brillouin zone is known, the paper details whether or not the 
Bloch functions of particular energy bands can be unitarily transformed into best localized Wannier 
functions symmetry-adapted to the space group G, to the magnetic group M, or to a subgroup of 
G or M. In this context, the paper considers usual as well as spin-dependent Wannier functions, 
the latter representing the most general definition of Wannier functions. The presented group 
theory is a review of the theory published by one of the authors in several former papers and is 
independent of any physical model of magnetism or superconductivity. However, it is suggested 
to interpret the special symmetry of the best localized Wannier functions in the framework of a 
nonadiabatic extension of the Heisenberg model, the nonadiabatic Heisenberg model. On the basis 
of the symmetry of the Wannier functions, this model of strongly correlated localized electrons makes 
clear predictions whether or not the system can possess superconducting or magnetic eigenstates. 
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I. INTRODUCTION 

The picture of strongly correlated localized or nearly- 
localized electrons is the base of a successful theoreti¬ 
cal description of both high-temperature superconductiv¬ 
ity and magnetism (see, e.g., [1-3] and citations given 
there). In almost all cases the appertaining localized elec¬ 
tron states are represented by atomic orbitals that define, 
for instance, partially filled s-, d -, or p- bands. 

Another option would be to represent the localized 
electron states by best localized and symmetry-adapted 
Wannier functions. In contrast to atomic functions, Wan¬ 
nier functions situated at adjacent atoms are orthogonal 
and, hence, electrons occupying (temporarily) adjacent 
localized states represented by Wannier functions com¬ 
ply with the Pauli principle. In addition, Wannier func¬ 
tion form a complete set of basis functions within the 
considered narrow, partially filled band. Consequently, 
Wannier functions contain all the physical information 
about this energy band. 

Wannier functions tend to be ignored by the theory 
of superconductivity and magnetism because we need a 
closed complex of energy bands [Definition 2.1] for the 
construction of best localized Wannier functions. Such 
closed complexes, however, do not exist in the band struc¬ 
tures of the metals where all the bands are connected to 
each other by band degeneracies. 

Fortunately, this problem can be solved in a natural 
way by constructing Wannier functions with the reduced 
symmetry of a magnetic group or by constructing spin- 
dependent Wannier functions as shall be details in the 
present paper. In both cases, interfering band degenera¬ 
cies are sometimes removed in the band structure with 
the reduced symmetry. 

Against the background of the described characteris¬ 


tics of the Wannier functions, our following two observa¬ 
tions should not be too surprising: 

(i) Materials possessing a magnetic structure with the 
magnetic group M also possess a closed, narrow 
and roughly half-filled complex of energy bands in 
their band structure whose Bloch functions can be 
unitarily transformed into best localized Wannier 
functions that are symmetry-adapted to the mag¬ 
netic group M. These energy bands form a “mag¬ 
netic band”, see Definition 6.2. 

(ii) Both normal and high-temperature superconduc¬ 
tors (and only superconductors) possess a closed, 
narrow and roughly half-filled complex of en¬ 
ergy bands in their band structure whose Bloch 
spinors can be unitarily transformed into best lo¬ 
calized spin-dependent Wannier functions that are 
symmetry-adapted to the (full) space group G of 
the material. These energy bands form a “super¬ 
conducting band”, see Definition 7.6. 

The first observation (i) was made at the band struc¬ 
tures of Cr [4], Fe [5], La 2 Cu 04 [6], YBa 2 Cu 3 C >6 [7], un¬ 
doped LaFeAsO [8], and BaFe 2 As 2 [9]; the second obser¬ 
vation (ii) at the band structures of numerous elemental 
superconductors [10] and of the (high-tenrperature) su¬ 
perconductors La 2 Cu 04 [6], YBa 2 Cu 3 C >7 [11], MgB 2 [11], 
and doped LaFeAsO [12]. It is particularly important 
that partly filled superconducting bands cannot be found 
in those elemental metals (such as Li, Na, K, Rb, Cs, Ca 
Cu, Ag, and Au) which do not become superconduct¬ 
ing [10]. An investigation into the band structures of 
the transition metals in terms of superconducting bands 
straightforwardly leads to the Matthias rule [13]. 

Though these two observations are clear, their theoret¬ 
ical interpretation is initially difficult. This is primarily 
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due to the fact that the models of localized electrons de¬ 
veloped so far, as, e.g., the familiar Hubbard model [14], 
are tailored to atomic orbitals that represent the localized 
states during an electronic hopping motion. Within mod¬ 
ern theoretical concepts, the Wannier functions often are 
nothing but a complete basis in the space spanned by the 
Bloch functions. Thus, their symmetry is often believed 
to do not tell anything about the physics of strongly cor¬ 
related electrons. 

In the light of this background, we suggest to interpret 
the special symmetries of best localized Wannier func¬ 
tions within the nonadiabatic Heisenberg model [11, 15]. 
This model of strongly correlated localized electrons 
starts in a consistent way from symmetry-adapted and 
best localized Wannier functions that represent the lo¬ 
calized electron states related to the hopping motion and 
defines the Hamiltonian H n of the related nonadiabatic 
system. On the basis of the symmetry of the Wannier 
functions, the nonadiabatic model makes clear predic¬ 
tions whether or not H n can possess superconducting or 
magnetic eigenstates [4, 9, 10, 16]. In this context, the 
nonadiabatic Heisenberg model no longer uses terms like 
s-, p-, or d-bands, but only speaks of superconducting or 
magnetic bands. 

In particularly interesting cases, the nonadiabatic 
Heisenberg model predicts that a small distortion of the 
lattice or a doping is required for the stability of the 
superconducting or magnetic state. Thus, in undoped 
LaFeAsO [8] and in BaFe 2 As 2 [9] the antiferromagnetic 
state must be stabilized by an experimentally well es¬ 
tablished distortion [17, 18], while in YBa2Cu306 [7] 
it is stable in the undistorted crystal. Superconduc¬ 
tivity in LaFeAsO [8] requires the experimentally con¬ 
firmed doping [18-21]. Also the superconducting state 
in LiFeAs [22] should be accompanied by a small dis¬ 
tortion of the lattice which, to our knowledge, is ex¬ 
perimentally not yet confirmed. Superconductivity in 
YBa 2 Cu 3 0 7 [11], MgB 2 [11] as well as in the transition 
elements [10] (such as in Nb [16]), on the other hand, 
does not require any distortion or doping. 

In the case of (conventional and high-T c [23]) supercon¬ 
ductivity, the nonadiabatic Heisenberg model provides a 
new mechanism of Cooper pair formation which may be 
described in terms of constraining forces [16] and spring- 
mounted Cooper pairs [24]. 

Any application of the nonadiabatic Heisenberg model 
starts with a determination of the symmetry of best local¬ 
ized (spin-dependent) Wannier functions related to the 
band structure of the material under consideration. In 
the following we shall summarize and update the group 
theory of Wannier functions as published in former pa¬ 
pers and give a detailed description how to determine 
the symmetry of best localized Wannier functions if they 
exist in the given band structure. Though we shall also 
define the two terms “magnetic” and “superconducting” 
bands which are related to the nonadiabatic Heisenberg 
model, the presented group theory is independent of any 
physical model of magnetism or superconductivity. 


II. USUAL (SPIN-INDEPENDENT) WANNIER 
FUNCTIONS 

II.1. Definition 

Consider a closed complex of p energy bands in the 
band structure of a metal or a semiconductor. 

Definition 2.1 (closed). A complex of energy bands is 
called closed if the bands are not connected by degenera¬ 
cies to bands not belonging to the complex. 

Definition 2.2 (closed band). In the following a closed 
complex of p energy bands is referred to as a single closed 
band consisting of p branches. 

The metals do not possess closed bands in their band 
structures. However, closed bands may arise after the ac¬ 
tivation of a perturbation reducing the symmetry in such 
a way that interfering band degeneracies are removed. 
Such a reduction of the symmetry may be caused by a 
magnetic structure or by a (slight) distortion of the crys¬ 
tal. 

Hence, we assume that the Hamiltonian H of a single 
electron in the considered material consists of a part Hg 
with the unperturbed space group G and a perturbation 
TLh with the space group H , 

H = 'H g +'H h , (2.1) 

where H is a subgroup of G, 

H C G. (2.2) 

In general, the considered closed energy band of p 
branches was not closed before the perturbation Hh was 
activated. 

Assume the Bloch functions <Pk,q(r) (labeled by the 
wave vector k and the branch index q) as the solutions 
of the Schrodinger equation of H to be completely calcu¬ 
lated in the first domain of the Brillouin zone. 

Definition 2.3 (first domain). Let be h the order of the 
point group Hq of H. Then the Brillouin zone is divided 
by the planes of symmetry into h domains. An arbitrary 
chosen domain we call the first domain. This first do¬ 
main shall comprise the bounding planes, lines and points 
of symmetry, too. 

As in Ref. [25], in the rest of the Brillouin zone the 
Bloch functions shall be determined by the equation 

Tak,q(r) = P({a\t a })ipk,q(r) for a <E H 0 , (2.3) 

where k lies in the first domain, and in the k space out¬ 
side the Brillouin zone by the equation 

Vk+K, q (r) = Tk, q (r). (2.4) 

K denotes a vector of the reciprocal lattice and H 0 stands 
for the point group of H. 
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Definition 2.4 (symmetry operators). P(a) denotes the 
symmetry operator assigned to the space group operation 
a = {a\t a } consisting of a point group operation a and 
the associated translation t a , acting on a wave function 
f(r) according to 

P(a)f(r) = /(a _1 r) = f(a~ 1 r - a~H a ). (2.5) 

The Bloch functions v 3 fc,g( T ’) of the closed band under 
observation can be unitarily transformed into Wannier 
functions 

BZ 

Wi(r-R- pi) = —j= V e~ zk{R+Pi) lp k4 (r) (2.6) 
v N k 

centered at the positions R + pi , where the functions 

Tk,i(r) = 9iq(k)<p k ,q{r) (2.7) 

9=1 

are “generalized” Bloch functions [25]. The sum in 
Eq. (2.6) runs over the N vectors k of the first Brillouin 
zone (BZ), the sum in Eq. (2.7) over the p branches of 
the considered band, R denote the vectors of the Bravais 
lattice, and the coefficients gi q (k) in Eq. (2.7) are the 
elements of an unitary matrix g (k), 

g _1 (fc) = g f (fc), (2.8) 

in order that the Wannier functions are orthonormal, 

J w*{r - R- Pi)wi'(r - R! - pr)dr = 5 RR '5n'- (2.9) 

Definition 2.5 (best localized). The Wannier functions 
are called best localized if the coefficients gi q {k ) may be 
chosen in such a way that the generalized Bloch functions 
<Pk,i(r) move - for fixed r - continuously through to whole 
k space [25]. 

As it was already shown in Ref. [26], the Bloch func¬ 
tions (pk, q {r) as the eigenfunctions of the Hamiltonian TL 
may be chosen in such a way that they vary continuously 
as functions of k through the first domain and, in par¬ 
ticular, approach continuously the boundaries of the first 
domain. From Eqs. (2.3) and (2.4), however, we cannot 
conclude that they also cross continuously the bound¬ 
aries of the domains within the Brillouin zone or at the 
surface of the Brillouin zone. Fortunately, this problem 
is solvable by group-theoretical methods [25, 27]. Theo¬ 
rem 4.1 shall define the condition for best localized and 
symmetry-adapted [Definition (2.7)] Wannier functions. 

II.2. Symmetry-adapted Wannier functions 

In Ref. [25] we demanded that symmetry-adapted 
Wannier functions satisfy the equation 

v 

w i (a~ 1 (r-R-p i )) = ^ Djj ( a)wj(r-R-pf) (2.10) 

3 =1 


for the elements a of the point group Hq of H, where the 
Dji(a) are the elements of the matrices 

D(a) = [A»] (2.11) 

forming a representation D of Hq which in most cases is 
reducible. [It should be noted that the sum in Eq. (2.10) 
runs over Wj(r — R— pi) and not over Wj(r — R — Pj)-\ 
Eq. (2.10) defines the symmetry of Wannier functions 
in general terms, particularly they may be centered at a 
variety of positions pi being different from the positions 
of the atoms. However, in the context of superconducting 
and magnetic bands we may restrict ourselves to Wannier 
functions centered at the positions of the atoms. 

Thus, we assume 

(i) that the positions pi of the Wannier functions in 
Eq. (2.6) are the positions of atoms, 

(ii) that only atoms of the same sort are considered 
(although, of course, other atoms may exist), and 

(iii) that there is one Wannier function at each atom. 
Under these assumptions [15], 

— the Wannier functions may be labeled by the posi¬ 
tions of the atoms, 

w T {r) = Wi(r - R-p^, (2.12) 

where 

T = R + Pi , (2.13) 

— the matrix representatives D(a) of the representa¬ 
tion D in Eq. (2.10) have one non-vanishing ele¬ 
ment Dij(a) with 

|A:»| = 1 (2.14) 

in each row and each column, and 

— Eq. (2.10) may be written in the considerably sim¬ 
plified form 

P(a)wr(r) = Dji(a)wT'(r) for a € H (2.15) 
where 

T' = aT + t a (2.16) 

and the subscripts i and j denote the number of 
the atoms at position T and T' , respectively. 

Definition 2.6 (number of the atom). The subscript i 
of the vector pi in Eq. (2.13) defines the number of the 
atom at position T. 

Definition 2.7 (symmetry-adapted). We call the Wan¬ 
nier functions symmetry-adapted to H if they satisfy 
Eq. (2.15). 
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Theorem 2.1. The third assumption (Hi) shows imme¬ 
diately that the number p of the branches of the band 
under observation equals the number of the considered 
atoms in the unit cell. 


Eqs. (2.15) and (2.16) define the non-vanishing ele¬ 
ments and, hence, we may write Eq. (2.14) more pre¬ 
cisely, 




1 if a pi +t a = pj + R 
0 else, 


(2.17) 


where {a\t a } £ H and R still denotes a lattice vector. 

Definition 2.8 (the representation defining the Wannier 
functions). In what follows, the representation D of Hq 
with the matrix representatives 


D(a) = [Dij(a)\ 

defined by Eq. (2.15) shall be shortly referred to as “the 
representation defining the Wannier functions” and its 
matrix representatives D(a) to as “the matrices defining 
the Wannier functions”. 

Definition 2.9 (unitary generalized permutation matri¬ 
ces). Since the matrices D(a) defining the Wannier func¬ 
tions have one non-vanishing element obeying Eq. (2.17) 
in each row and each column, they are so-called unitary 
generalized permutation matrices. 


III. DETERMINATION OF THE 
REPRESENTATIONS D DEFINING THE 
WANNIER FUNCTIONS 


In the following Sec IV we shall give a simple con¬ 
dition [Theorem 4.1] for best localized and symmetry- 
adapted Wannier functions yielding the representations 
of the Bloch functions at all the points k of symmetry 
in the Brillouin zone. However, in Theorem 4.1 the rep¬ 
resentations D defining the Wannier functions must be 
known. Hence, first of all we have to determine in this 
section all the possible representations that may define 
the Wannier functions. In this context we assume first 
that all the atoms are connected by symmetry. This 
restricting assumption shall be abandoned not until in 
Sec. III.4. 

Definition 3.1 (connected by symmetry). Two atoms at 
positions pi and pj are connected by symmetry if there 
exists at least one element a = {a\t a } in the space group 
H satisfying the equation 


api + t a — pj + R, (3.1) 


III.1. General properties of the representatives 
D(a) of D 


First consider the diagonal elements 

di(a) = Du(a) (3.2) 

of the matrices D(a) defining the Wannier functions. 
From Eq. (2.17) we obtain 

17 / \i / 1 if Oipi -\- t a — pi -f- R /q q\ 
= { 0 else (3 ' 3) 

where R denotes a lattice vector. This equation demon¬ 
strates that the matrix D(a) has non-vanishing diagonal 
elements di(a) if the space group operation a = {a\t a } 
leaves invariant the position pi of the ith atom. These 
space group operations form a group, namely the group 
G Pi of the position pi. 

Definition 3.2 (group of position). The group G Pi of 
the position pi is defined by 

a £ G Pi if a £ H and api + t a = pi + R. (3.4) 

Go Pi denotes the point group of G Pi . 

Hence, the non-vanishing diagonal elements di(a) of 
the matrices D(a) form a one-dimensional representation 
di of the point group Go Pi of G Pi . The Wannier functions 
transform according to 

P(a)w T {r) = di(a)w T +n.(r) for a £ G 0pi (3.5) 

[cf. Eq. (2.15)] by application of a space group opera¬ 
tor P(a) (where R still denotes a vector of the Bravais 
lattice). From Eq. (2.10) we may derive the equivalent 
equation 

Wi(a~ l (r-R-pi)) = di(a)wi(r - R- p^ for a £ G 0pi 

(3-6) 

or, after shifting the origin of the coordinate system into 
the center of the function Wi(r — R— pi), 

r' = r - R - pi, 

we receive an equation 

w i (a _ 1 r') = di(a)wi(r') for a £ G 0pi (3.7) 

emphasizing the point-group symmetry of the Wannier 
function at position R + pi. 

In constructing the representation D defining the 
Wannier functions we cannot arbitrarily chose the one¬ 
dimensional representations di of Go Pi because they must 
be chosen in such a way that the matrix representatives 
D (a) form a representation of the point group Hq, i.e., 
they must obey the multiplication rule 

D(a/3) = D(a)D(/3) (3.8) 

for all the elements a and /3 in Hq. 


where R is a lattice vector. 
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In what follows we assume that all the groups G Pi are 
normal subgroups of H. In fact, in all the crystal struc¬ 
tures we examined in the past, G Pi was a normal sub¬ 
group, be it because it was a subgroup of index 2 or be 
it because it was the intersection of two subgroups of in¬ 
dex 2. Both cases are sufficient for a normal subgroup. 
We believe that in all physically relevant crystal struc¬ 
tures G Pi is a normal subgroup of H. If not, the present 
formalism must be extended for these structures. 

When the groups G Pi are normal subgroups of H , each 
of the groups G Pt contains only complete classes of H , 

b~ 1 ab £ G Pi if a £ G Pi and b £ H. (3.9) 

We now show that, as a consequence, all the groups G Pi 
contain the same space group operations. 

Let be b = {f3\tp} a space group operation of H moving 
pi into pj, 

flpi +tp = pj + R , 

then 

c = b~ 1 ab (3.10) 

is an element of G Pi if a £ G Pj . Eq. (3.10) even yields all 
the elements c of G Pi when a runs throw all the elements 
of G Pj because we may write Eq. (3.10) in the form 

bob- 1 = a 

showing that we may determine from any element c £ 
G Pi an element a £ G Pj . 

On the other hand, Eq. (3.9) shows that c is an element 
of G p . , too. When a runs through all the elements of 
G Pj , then also c runs through all the elements of G Pj . 
Consequently, all the groups G Pi as well as all the related 
point groups Go Pi contain the same elements. 

Thus, we may omit the index i and define 

Definition 3.3 (group of position). The group G p and 
the related point group Gq p of the positions of the atoms 
is defined by 

G p = G Pi (3.11) 

and 

G 0p = G 0pi , (3.12) 

respectively, where G Pi and Go Pi are given by Defini¬ 
tion 3.2. 


III.2. Necessary condition for of the 
representatives D(a) of D 

The one-dimensional representations di of Go P must 
be chosen in such a way that the matrices D(a) defining 
the Wannier functions form a representation D of the 
complete point group Hq. A necessary condition is given 
by the evident Theorem 3.1. 


Theorem 3.1. If the matrices D(a) cannot be com¬ 
pletely reduced into the irreducible representations of Ho, 
then they do not form a representation of the point group 
Ho- 

This theorem is necessary, but not sufficient: even if 
the matrices D(a) can be completely reduced into the 
irreducible representations of Hq then they need not form 
a representation of the point group H 0 [28]. The complete 
decomposition of a reducible representation is described, 
e.g., in Refs. [28] and [29], in particular, see Eq. (1.3.18) 
of Ref. [29]. Theorem 3.1 leads to three important cases: 

— Case (i): If all the representations d, are subduced 
from one-dimensional representations of Hq, then 
all the representations d, are equal, 

di = d for all the positions pi. (3.13) 

The representation d may be equal to any one¬ 
dimensional representation of Gq p subduced from 
a one-dimensional representation of Hq. 

— Case (ii): If all the representations di are subduced 
from two-dimensional representations of Ho, then 
one half of the representations di is equal to d^ 
and the other half is equal to dg, 

di = dA for one half of the positions pi 

di = dg for the remaining positions pi, - 

where dA and dg are subduced from the same two- 
dimensional representation of Hq. In special cases, 
the two representations d^ and dg may be equal, 
see below. 

— Case (iii): “Mixed” representations D consisting 
of both representations d; subduced from one- and 
two-dimensional representations of Hq do not exist. 

A further case that the representations di are subduced 
from f/iree-dimensional representations of Hq may occur 
in crystals of high symmetry but is not considered in this 
paper. 

These results (i) - (iii) follow from the very fact that 
Eq. (2.15) describes an interchange of the Wannier func¬ 
tions at different positions pi. Such an interchange, how¬ 
ever, does not alter the symmetry of the Wannier func¬ 
tions. 


III.3. Sufficient condition for of the representatives 
D(a) of D 

For a £ Go P the matrices D(a) defining the Wan¬ 
nier functions are diagonal, while the remaining matrices 
D( a) [for a £ H 0 — Gq p ] do not possess any diagonal 
element. Theorem 3.1 only gives information about the 
diagonal matrices D(cc). Hence, this theorem indeed can¬ 
not be sufficient because we do not know whether or not 
the remaining matrices obey the multiplication rule (3.8). 



6 


In this section we assume that the matrices D(a) al¬ 
ready satisfy Theorem 3.1 and examine the conditions 
under which they actually form a (generally reducible) 
representation of Hq. In doing so, we consider separately 
the two cases (i) and (ii) of the preceding Sec. III.2. 


III. 3.1. Case (i) of Sec. III. 2 

No further problems arises when case (i) of Sec. III.2 is 
realized. In this case, Theorem 3.1 is necessary and suffi¬ 
cient. To justify this assertion, we write down explicitly 
the non-diagonal elements of the matrices D(a). 

Let be 6 any one-dimensional representation of Ho sub- 
ducing the representation d in Eq. (3.13). If we put all 
the non-vanishing elements of the matrices D(a) equal 
to the elements 5(a) of 5, 

Dji(a) = { S Q {a) «£* + t « = ?i+ R (3.15) 

then we receive matrices D(a) evidently multiplying as 
the elements of the representation S and, consequently, 
obeying the multiplication rule in Eq. (3.8). 


III. 3.2. Case (ii) of Sec. III. 2 

The situation is a little more complicated when case 
(ii) of Sec. III.2 is realized. Now, the representations eLi 
and ds in Eq. (3.14) may be distributed across the po¬ 
sitions pi in such a way that the matrices D(a) form a 
representation of Hq or do not. Though we always find a 
special distribution of the dA and ds yielding matrices 
D (a) actually forming a representation of Ho, we have to 
rule out those distributions not leading to a representa¬ 
tion of Ho, because in the following [in Eqs. (4.1), (6.10), 
and (7.42)] we need the matrices D(a) explicitly. 

Let be A [with the matrix representatives A (a)] a 
two-dimensional representation of Ho subducing the two 
representations dA and ds of Gq p ■ The matrix represen¬ 
tatives A (a) may be determined, e.g., from Table 5.1 of 
Ref. [29], 

As a first step, A must be unitarily transformed (by 
a matrix Q) in such a way that the matrices A (a) are 
diagonal for a € Gq p , 

A(a) = Q _1 A(a)Q = diagonal for a € Go P ■ (3.16) 

Now consider a certain distribution of the representa¬ 
tions dA and ds across the positions Pi- Then we may 
determine the elements of the matrices D(a), if they ex¬ 


ist, be means of the formula 
if apt + t a = pj + R 

Ai 2 (a) if dj = dA and d, = dg, 

A 2 i(a) if dj = ds and d, = dA, 

= < _ 

An(a) if dj = dA and d, = d^, 

A 2 2 (a) if dj = ds and d, = ds, 

else 

Dji(a) = 0, 

(3.17) 

where the A ij(a) denote the elements of A (a). 

It turns out that in each case the matrices determined 
by Eq. (3.17) satisfy the multiplication rule in Eq. (3.8) if 
Eq. (3.17) produces for each space group operation a € H 
an unitary generalized permutation matrix D(a). This 
may be understood because Eq. (3.17) defines the com¬ 
plex numbers Dji(a) in such a way that the Wannier 
functions transform in Eq. (2.15) in an unequivocal man¬ 
ner like the basis functions for A. With “like” the basis 
functions we want to express that by application of any 
space group operator P({a\t a }) they are multiplied in 
Eq. (2.15) by the same complex number A ij(a) as the 
basis functions for A. The Wannier functions would in¬ 
deed be basis functions for A if they would not be moved 
from one position pi to another by some space group op¬ 
erations. Hence, we may expect that the matrices D(a) 
satisfy the multiplication rule in Eq. (3.8) just as the 
matrices A (a) do. Nevertheless, the multiplication rule 
should be verified numerically in any case. 

When using this Eq. (3.17) a little complication arises 
if the group of position Go P contains so few elements 
that the two one-dimensional representations d\ and ds 
subduced from A are equal. Thus, in this case we have 
no problem with the distribution of d\ and ds across 
the positions pi. Theorem 3.1 is necessary and sufficient 
and we may directly solve Eq. (4.1) of Theorem 4.1. 

However, when in Sec. VI or in Sec. VII.3 we will con¬ 
sider magnetic groups, we need all the representatives 
D (a) of the representation D explicitly. Fortunately, also 
when the representations dA and ds are equal, Eq. (3.17) 
is applicable: in this case their exists at least one diago¬ 
nal matrix representative A ( 7 ) of A with vanishing trace 
and 7 ^ Gq p - We may define pairs 

( Pa,Pb), (Pc, Pd), (3.18) 

of positions pi where the positions in each pair are con¬ 
nected by the space group operation { 7 |i 7 }- In the sim¬ 
plest case, we receive two pairs. Then in Eq. (3.17) we 
may identify the two representations at p a and Pb by 
dA and the representations at the other two positions 
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p c and pd by dg- If we find four pairs of positions, we 
may look for a second matrix representative A(y / ) in A 
with vanishing trace and 7 ' ^ Go p . Then we may repeat 
the above procedure and receive again four pairs of posi¬ 
tions. Now we associate the two representations dA and 
cLb to the positions pi under the provision that always 
positions of the same pair are associated with the same 
representation dn or d^. 

Finally, it should be mentioned that the elements of 
the non-diagonal matrices D(a) are not fully fixed (as 
already remarked in Ref. [27]): In Eq. (3.15) we may use 
the elements <5(a) of any one-dimensional representation 
S subducing the representation d. We receive in each 
case the same diagonal, but different non-diagonal ma¬ 
trices nevertheless satisfying the multiplication rule (3.8). 
Analogously, in Eq. (3.17) we may determine the matrices 
D(cc) by means of any two-dimensional representation A 
subducing dn and d^. 

In the following Theorem 3.2 we summarize our results 
in the present Sec. III.3. 


Theorem 3.2. The Wannier function wi{r R pi) 
at the position pi is basis function for a one-dimensional 
representation di of the “point group of position” Gq p C 
Hq [Definition 3.3], cf. Eq. (3.7). The representations di 
fix the (generally reducible) representation D of Hq defin¬ 
ing the Wannier functions [Definition 2.8]. The matrix 
representatives D(a) of D are unitary generalized per¬ 
mutation matrices. We distinguish between two cases (i) 
and (ii). 

Case (i): If the representations di are subduced from 
one-dimensional representations of the point group Hq, 
then all the Wannier functions of the band under obser¬ 
vation are basis functions for the same representation d 
which may be any one-dimensional representation of Gq p 
subduced from a one-dimensional representation of H 0 . 
The representation D exists always, its matrix represen¬ 
tatives D(a) may be calculated by Eq. (3.15). 

Case (ii): If the representations di are subduced from 
two-dimensional representations of the point group H 0 , 
then the Wannier functions are basis functions for the 
two one-dimensional representations dA and dg of Gq p 
subduced from the same two-dimensional representation 
of Hq ■ One half of the Wannier functions is basis func¬ 
tion for dA and the other half for dg. In special cases, 
the representations dA and ds may be equal, see above. 
The representation D exists for a given distribution of 
the representations d \ and ds across the positions pi if 
Eq. (3.17) yields unitary generalized permutation matri¬ 
ces D(a) satisfying the multiplication rule in Eq. (3.8). 

A third case with representations di subduced from one¬ 
dimensional as well as from two-dimensional representa¬ 
tions of Hq does not exist. 


III.4. Not all the atoms are connected by 
symmetry 


If not all the atoms at the positions pi are connected by 
symmetry [Definition 3.1], the representation D defining 
the Wannier functions consists of representatives D(a) 
which may be written in block-diagonal form 


D(a) = 


block 1 


block 2 • • • 


(3.19) 


where each block comprises the matrix elements Dij{a) 
belonging to positions connected by symmetry. Other¬ 
wise, when the matrices D(a) would not possess block- 
diagonal form, Eq. (2.10) would falsely connect atomic 
positions that are not at all connected by symmetry. As 
a consequence of the block-diagonal form, the representa¬ 
tion D is the direct sum over representations D q related 
to the individual blocks, 

D = D 1 © D 2 © ... 

= Y^D q . (3.20) 

We may summarize as follows. 

Theorem 3.3. Each block D q in Eq. (3.20) forms its 
own representation of Hq and, hence, must comply sepa¬ 
rately and independently with the criteria given in Theo¬ 
rem 3.2. 


The groups of position G p belonging to different blocks 
may (but need not) be different. However, we assume 
that the sum in Eq. (3.20) consists only of blocks with 
coinciding groups of position. If this is not true in special 
cases, the number p of the atoms in Eq. (2.7) must be 
reduced until the groups of position coincide in the sum 
in Eq. (3.20). Briefly speaking, in such a (probably rare) 
case atoms of the same sort must be treated like different 
atoms. 


IV. CONDITION FOR BEST LOCALIZED 
SYMMETRY-ADAPTED WANNIER FUNCTIONS 

Remember that we consider a closed energy band of 
p branches and let be given a representation D defining 
the Wannier functions which was determined according 
to Theorems 3.2 and 3.3. Then we may give a simple 
condition for best localized symmetry-adapted Wannier 
functions based on the theory of Wannier functions pub¬ 
lished in Refs. [25] and [27]. 

Theorem 4.1. Let be k a point of symmetry in the first 
domain of the Brillouin zone for the considered material 



and let be Hk C H the little group of k in Herrings sense. 
That means, H & is the FINITE group denoted in Ref. [29] 
by H G k (and listed for all the space groups in Table 5.7 
ibidem). Furthermore, let be Dk the p-dimensional rep¬ 
resentation of Hk whose basis functions are the p Bloch 
functions pk,q(r) with wave vector k, and X k{o.) (with 
a € Hk) the character of Dk ■ Dk either is irreducible 
or the direct sum over small irreducible representations 

of H k - 

We may choose the coefficients gi q (k ) in Eq. (2.7) in 
such a way that the Wannier functions are best local¬ 
ized [Definition 2.1] and symmetry-adapted to H [Defi¬ 
nition 2.7] if the character Xk{a) of Dk satisfies at each 
point k of symmetry in the first domain of the Brillouin 
zone the equation 

Xk (a) = e~ iak t ° for a € H k , 

1=1 

(4.1) 

where a = {a\t a } and 

n i (a) = |o i(a) if dse a€G ° P (4-2) 

The complex numbers di(a) stand for the elements of 
the one-dimensional representations d, of Gq p fixing the 
given p-dimensional representation D defining the Wan¬ 
ner functions. 

Definition 4.1 (point of symmetry). The term “point 
of symmetry” we use as defined in Ref. [29]: k is a point 
of symmetry if there exists a neighborhood of k in which 
no point except k has the symmetry group Hk ■ 

Thus, a point k of symmetry has a higher symmetry 
than all surrounding points. 

We add a few comments on Theorem 4.1. 

— In Eq. (4.2) we write n,;(a) rather than rii(a) be¬ 
cause the group Go p depends on a = {a\t a }. 

— The representation D defining the Wannier func¬ 
tions is equivalent to the representation Do, i.e., to 
the representation Dk for k = 0, see Eq. (5.10). 

— In the majority of cases all the representations d, 
in Eq. (4.2) are equal. The only exceptions arises 
when 

(i) not all the positions pi are connected by sym¬ 
metry or 

(ii) the one-dinrensional representations d, of Go p 
are subduced from a higher-dimensional rep¬ 
resentation of Hq. 

— A basic form of Theorem 4.1 was published first 
in Eq. (23) of Ref. [ 6 ] and used in several former 
papers. Eq. (23) of Ref. [ 6 ] yields the same results 
as Theorem 4.1 

(i) if all the p; are connected by symmetry and 


(ii) if all the representations d t of Gq p are sub¬ 
duced from one-dimensional representations of 

H 0 . 

These two conditions were satisfied in our former 
papers. 

— The irreducible representations of the Bloch func¬ 
tions of the considered band at the points k of sym¬ 
metry may be determined from the representations 
D k as follows: 

Theorem 4.2. Let Hk possess r irreducible represen¬ 
tations with the characters X k,m(o) (1 < m < r) and 
assume that Dk contains the mth irreducible representa¬ 
tion, say, c m times. Then the numbers c m may be calcu¬ 
lated by means of Eq. (1.3.18) of Ref. [29], 

1 H h 

Cm — | TT | Xk,m{o)Xk( a )i ( 4 - 3 ) 

\ Hk \ V 

where X k (a) denotes the character of Dk as determined 
by Eq. (f.l) and the sum runs over the \Hk\ elements 
a of Hk- Remember [Theorem f.l] that Hk is a finite 
group. 

V. PROOF OF THEOREM 4.1 

The existence of best localized symmetry-adapted 
Wannier functions is defined in Satz 4 of Ref. [25]: such 
Wannier functions exist in a given closed energy band of 
p branches if Eqs. (4.28) and (4.17) of Ref. [25] are satis¬ 
fied. We show in this section that the fundamental Theo¬ 
rem 4.1 complies with these two equations if the Wannier 
functions meet the assumptions (i) - (iii) in Sec. II.2. 

V.l. Equation (4.28) of Ref. [25] 

As a first step consider Eq. (4.28) of Ref. [25] stat¬ 
ing that best localized and symmetry-adapted Wannier 
functions may exist only if two representations Dy 1 , R and 
D k > are equivalent, 

Dk equivalent to Dk, (5-1) 

where we have abbreviated k'j, R by k denoting a point 
of symmetry lying in the first domain of the Brillouin 
zone. Consequently, our first task will be to determine 
the character of Dk as well as of Dk, 

The representation Dk as defined in Theorem 4.1 is the 
direct sum of the representations of the Bloch functions 
of the considered band at point k. The character Xk{a) 
of the representation Dk is simply given by 

Xfc(a) = trace D fe (a) (5.2) 

where the matrices D*. (a) are the matrix representatives 
of Dk- 
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The matrix representatives D*,(a) of D k are defined in 
Eq. (4.26) of Ref. [25], 


in Eq. (5.5) is a matrix representative D(a) of the rep¬ 
resentation D defining the Wannier functions, 


D fc (a) = S* (K a )D 0 (a)e~ iakta (5.3) 

where 

K a =k- ak (5.4) 

is a vector of the reciprocal lattice. Again we have 
abbreviated k' SR by k denoting a point of symmetry. 
The matrices S (K) as defined in Eq. (4.13) of Ref. [25] 
are responsible for a continuous transition of the gen¬ 
eralized Bloch functions between neighboring Brillouin 
zones. The matrices Do (a) are the matrix representa¬ 
tives of the representation Dk for k = 0 as defined in 
Theorem 4.1. Do is the direct sum of the irreducible 
representations of the Bloch functions of the considered 
band at point T. 

The traces of the matrices Dfc(a) can be determined 
by transforming Eq. (5.3) with the complex conjugate of 
the matrix M defined by Eq. (2.1) of Ref. [27], 

M*Dfe(a)M*~ 1 = M*S*(ZT a )M* -1 x 

M*D 0 (a)M* -1 x (5.5) 

^—ioik-tcx 

where a = {a\t a } still denotes an element of the space 
group H. By definition, the matrix M diagonalizes the 
matrices S (K) which is possible since all the S (K) com¬ 
mute. Thus, the first factor M*S*(K’ a )M* _1 in Eq. (5.5) 
is the diagonal matrix 

S*(K a ) =e~ lK ^, (5.6) 

where, according to Eq. (2.7) of Ref. [27], also T is a 

diagonal matrix with 

Ta = Pi . (5.7) 

Hence, the first factor in Eq. (5.5) may be written as 
M*S*{K a )M*~ 1 = S*(K a ) = 

0 ° \ 

0 0 

e -ip2-(k-ak) q 

0 e -ip!-(k-ak) j 


o-ip^-(k-ak) 


(5.8) 

Definition 5.1 (horizontal bar). In line with Ref. [27], 
we denote matrices transformed with M (or by a 
horizontal bar to indicate that these matrices belong to 
the diagonal matrices S {K). 

As shown in Ref. [27] (see Eqs. (2.18) and (2.19) of 
Ref. [27]), the second factor 

D 0 (a) = M*D 0 (a)M * -1 (5.9) 


D 0 (a) = D(a). (5.10) 

Thus, the matrices 

D fc (o) = M*Dfc(a)M * -1 

= S*(K a )D(a)e- iQfc ‘ t “ (5.11) 

are the matrix representatives of a representation Dk 
equivalent to Dk- 

The character of Dk may be easily determined: The 
diagonal elements di(a) of the matrices D(a) are fixed by 
Theorems 3.2 and 3.3. Since the matrix S (K a ) is diag¬ 
onal, the diagonal elements di(a) of the matrices Dfc(a) 
may be written as 

di(a) = e- iak ta di(a)e- ipi ^ k - ak ^ for a G H k , (5.12) 

where still a = {n|t Q }. The diagonal elements di(a ) 
of the matrices D(a) vanish if a ^ Gq p , see Eq. (3.3). 
Hence, the term on the right-hand side of Eq. (4.1) is 
the sum over the diagonal elements di(a), i.e., it is the 

trace of the matrices D fc (a). Consequently, if Eq. (4.1) 
is satisfied then condition (5.1) is true. 

Strictly speaking, in Ref. [25] we have proven that the 
condition (5.1) must be satisfied for the points of symme¬ 
try lying in the first domain on the surface of the Brillouin 
zone. Eq. (4.1) demands that in addition the represen¬ 
tation Dq is equivalent to the representation D which is 
evidently true, see Eq. (5.10). 


V.2. Equation (4.17) of Ref. [25] 


As a second step we show that Eq. (4.17) of Ref. [25] 
does not reduce the validity of Theorem 4.1 but this 
equation is satisfied whenever the assumptions (i) - (iii) 
in Sec. II.2 are valid. Taking the complex conjugate of 
Eq. (4.17) of Ref. [25] and transforming this equation 
with the matrix M* already used in Eq. (5.5), we receive 
the equation 

S *{aK) = D(a)S\K)D-\a)e- iaKt «, (5.13) 


cf. Eqs. (5.8) and (5.10), which must be satisfied for all 
a = {a\t a } G H and all the vectors K of the reciprocal 
lattice. 

Just as the matrix 


S *(K) 


0 —ip^K 

0 

0 

0 


0 0 \ 

0 0 

e -iP2--FC 0 

0 e -ipi K j 


(5.14) 


also the matrix D(a)S*(ff)D 1 (cc) in Eq. (5.13) is di¬ 
agonal with the same diagonal elements which, however, 
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may stand in a new order. In fact, if D.j, (a) ^ 0, the ele¬ 
ment e~ lpi ' K of S *{K) at position i stands at position j 
in the matrix D (a )S (/f)D _1 (a). Thus, from Eq. (5.13) 
we receive the p equations 

e -icxK p 3 = e -iK Pi . e -i*K-t a if D ..( a ) ± 0, (5.15) 

yielding p equations for the positions p,, 

Pj = api+t a + Rj if Dji(a) ^ 0, (5.16) 

where Rj is a lattice vector which may be different in 
each equation. In fact, these last p equations (5.16) are 
satisfied, see Eq. (2.17). 

VI. MAGNETIC GROUPS 

Assume a magnetic structure to be given in the con¬ 
sidered material and let be 

M = H + K{^\t}H (6.1) 

the magnetic group of this magnetic structure, where 

{tI t } e G (6.2) 

and K denotes the operator of time inversion acting on 
a function f(r) of position according to 

Kf(r) = f*{r). (6.3) 

We demand that the equation 

Kwi( 7 -1 (r - R - pi)) = ^2 NjiWjir - R - pf) (6.4) 

j=i 

is satisfied in addition to Eq. (2.10), where the matrix 
N = [fV,;j] is the representative of the anti-unitary sym¬ 
metry operation A'y in the co-representation of the point 
group 

M 0 =H 0 + K'yHo (6.5) 

of M derived from [29] the representation D of Hq defin¬ 
ing the Wannier functions. 

Still we assume that there is exactly one Wannier func¬ 
tion at each position p, , i.e., the three assumptions (i) - 
(iii) of Sec. II.2 remain valid. Thus [15], Eq. (6.4) may 
be written in the more compact form 

KP({~/\T})w T (r) = NjiW T >{r) ( 6 . 6 ) 

with 

T' = 7 T + t (6.7) 

and the subscripts i and j denote the number of the 
atoms at position T and T' , respectively, see Defini¬ 
tion 2.6. 


Definition 6.1 (symmetry-adapted to a magnetic 
group). We call the Wannier functions symmetry- 
adapted, to the magnetic group M if, in addition to 
Eq. (2.15), Eq. (6.6) is satisfied. 

Again (cf. Sec. II.2), Eq. ( 6 . 6 ) defines the non¬ 
vanishing elements of the Matrix N. Hence, also N has 
one non-vanishing element in each row and each column, 

™={;ir +T=ft+K («> 

As already expressed by Eq. (6.4), we only consider 
bands of p branches which are not connected to other 
bands also after the introduction of the new anti-unitary 
operation A'{ 7 |t}. That means that the considered band 
consists of p branches as well after as before the intro¬ 
duction of A'{ 7 |t}. Hence, the matrix N must satisfy 
the equations 

NN* = D( 7 2 ) (6.9) 

and 

D(a) = ND*( 7 " 1 a 7 )N “ 1 for a <E H 0 , (6.10) 

see Eq. (7.3.45) of Ref. [29]. Still the matrices D(a) 
are the representatives of the representation D of Hq 
defining the Wannier functions. In Ref. [29] Eq. (7.3.45) 
was established for irreducible representations. However, 
this prove in Sec. 7.3 ibidem shows that Eq. (7.3.45) holds 
for reducible representations, too, if Eq. (6.9) is satisfied. 

Assume Theorem 4.1 to be satisfied in the consid¬ 
ered energy band and remember that then the coeffi¬ 
cients gi q (k ) in Eq. (2.7) can be chosen in such a way 
that the Wannier functions of this band are best local¬ 
ized and symmetry-adapted to H. In Ref. [30] we have 
shown that the Wannier functions may even be chosen 
symmetry-adapted to the magnetic group M if Eq. (7.1) 
of Ref. [30], 

S(— 7 EO = DS(A^ 7 )S*(A:)DS- 1 (^ 7 )e- i ^", (6.11) 

is valid for each vector K of the reciprocal lattice (which 
should not be confused with the operator K of time in¬ 
version). The matrix S(AT) is defined in Eq. (4.13) of 
Ref. [25] and the matrix Do(A'q) is the representative of 
the symmetry operation K'y in the co-representation of 
Mo derived from the representation Do, i.e., from the 
representation Dy, for fc = 0 as introduced in Theo¬ 
rem 4.1. 

Transforming Eq. (6.11) with the matrix M* already 
used in Eq. (5.5) and using 

N = M*Do(A' 7 )M - 1 (Eq. (11.29) of Ref. [30]) 
S (K) = MS(K)M ' 1 = diagonal, Eq. (5.14) 

S*( 7 AT) = S(- 7 AT) (see Eq. (5.14)) 

( 6 . 12 ) 

we receive an equation 

S *( 7 AT) = N*S*(K)N*- 1 e- i7 - f<: ' T 


(6.13) 
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identical to Eq. (5.13) when we replace the space group 
operation {a|t Q } by {qlr} and D(a) by N*. In Sec. V.2 
we have shown that Eq. (5.13) is satisfied if the matri¬ 
ces D(a) follow Eq. (2.17). In the same way, Eq. (6.13) 
is true if the elements of N (as well as of N*) obey 
Eq. ( 6 . 8 ). Thus, Eqs. ( 6 . 8 ), (6.9), and (6.10) are the 
only additional conditions for the existence of best local¬ 
ized Wannier functions which are symmetry-adapted to 
the magnetic group M. 

We summarize the results of the present Sec. VI in 

Theorem 6.1. The coefficients gi q {k) in Eqs. (2.7) may 
be chosen in such a way that the Wannier functions 
are best localized [Definition 2.5] and even symmetry- 
adapted to the magnetic group M in Eq. (6.1) [Defini¬ 
tion 6.1] if, according to Theorem f.l, they may be cho¬ 
sen symmetry-adapted to H and if, in addition, there ex¬ 
ists a fi-dimensional matrix N satisfying Eqs. (6.8), (6.9) 
and (6.10). 

The representation D in Eqs. (6.9) and (6.10) is the 
representation defining the Wannier functions as used in 
Theorem 4.1. 

In most cases, we may put the non-vanishing elements 
ofN equal to 1. 

Definition 6.2 (magnetic band). If, according to Theo¬ 
rem 6.1, the unitary transformation in Eq. (2.6) may be 
chosen in such a way that the Wannier functions are best 
localized and symmetry-adapted to the magnetic group M 
in Eq. (6.1), we call the band under consideration [as de¬ 
fined by the representations Dk in Eq. (4-1)] a, “magnetic 
band related to the magnetic group M”. 

Within the nonadiabatic Heisenberg model, the exis¬ 
tence of a narrow, roughly half-filled magnetic band in 
the band structure of a material is a precondition for the 
stability of a magnetic structure with the magnetic group 
M in this material. However, the magnetic group M 
must be “allowed” in order that the time-inversion sym¬ 
metry does not interfere with the stability of the magnetic 
state [9]. 

VII. SPIN-DEPENDENT WANNIER 
FUNCTIONS 

VII.1. Definition 

Assume the Hamiltonian TL of a single electron in the 
considered material to be given and assume H to con¬ 
sist of a spin-independent part Hi and a spin-dependent 
perturbation H s , 

H = Hi+H s . (7.1) 

Further assume the Bloch spinors '<pk,q,s{r, t) as the exact 
solutions of the Schrodinger equation 

7~L'^fk,q,s{j‘ = Ek,q,s'l]k,q,s(j' ,t) (7.2) 

to be completely determined in the first domain of the 
Brillouin zone. Just as the Bloch functions, they are 


labeled by the wave vector k and the branch index q. In 
addition, they depend on the spin coordinate t, = ±i and 
are labeled by the spin quantum number s = ± 5 . 

Consider again a closed energy band of /i branches 
which, in general, was not closed before the perturba¬ 
tion H s was activated. Now each branch is doubled, that 
means that it consists of two bands related to the two 
different spin directions. Just as in Sec. II. 1 we assume 
that the Bloch spinors i[)k,q,s{i‘H) are chosen in such a 
way that they vary continuously through the first do¬ 
main and approach continuously the boundaries of the 
first domain. In the rest of the Brillouin zone and in the 
remaining k space they shall be given again by Eqs. (2.3) 
and (2.4) [30] where, however, P(a) acts now on both r 
and t, see Eq. (7.16). 

We define “spin-dependent Wannier functions” by re¬ 
placing the Bloch functions <pk. q (r) i n Eq. (2.7) by linear 
combinations 

+ 5 

t)= E fms (q,k)i/j k , q ,s(r,t) (7.3) 
—;§ 

of the given Bloch spinors. Thus, Eq. (2.7) becomes 

*) = £ giq(k)lfi k ,q,m 0 r,t ) (7.4) 

9=1 

and, finally, the spin-dependent Wannier functions my be 
written as 

BZ 

Wi,m{r - R~ Pi,t) = -t=E ]e~ lk{R+Pi) Tk,i,m{r,t). 

^k 

(7.5) 

Also the spin-dependent Wannier functions depend on 
t and are labeled by a new quantum number rn = ±5 
which, in the framework of the nonadiabatic Heisen¬ 
berg model, is interpreted as the quantum number of the 
“crystal spin” [15, 31, 32]. The sum in Eq. (7.4) runs over 
the /i branches of the given closed energy band, where /1 
still is equal to the number of the considered atoms in 
the unit cell. 

The matrices 

g i k ) = [g iq {k )] (7.6) 

still are unitary [see Eq. (2.8)] and also the coefficients 
f ms {q,k) in Eq. (7.3) form for each k and q a two- 
dimensional matrix 

f {<l,k) = [fms(q,k)\ (7.7) 

which is unitary, 

f~ 1 (q,k) = f'(q,k), (7.8) 

in order that the spin-dependent Wannier functions are 
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orthonormal, 

+ 5 , 

y, / w* m (r - R- pi,t)wi’ tm >(r ~ R' pi>,t)dr 


Within the nonadiabatic Heisenberg model we strictly 
consider the limiting case of vanishing spin-orbit cou¬ 
pling, 

U s 0, (7.10) 

by approximating the Bloch spinors 4’k,q,s( r i t) by means 
of the spin-independent Bloch functions <fik,q(i')- 1 11 
this context, we should distinguish between two kinds 
of Bloch states ipk, q (i') in the considered closed band: 

(i) If £k, q {r) 

— was basis function for a non-degenerate rep¬ 
resentation already before the spin-dependent 
perturbation R s was activated, or 

— was basis function for a degenerate representa¬ 
tion before R s was activated, and this degen¬ 
eracy is not removed by H s [see Sec. VII.4.2], 

then we may approximate the Bloch spinors by 

fl>k,q,s{r,t) = u 8 (t)<pk,q{r) (7.11) 

where the functions u s {t) denote Pauli’s spin func¬ 
tions 

u s (t) = 5 st , (7.12) 

with the spin quantum number s = and the 
spin coordinate t = ±|. Eq. (7.11) applies to the 
vast majority of points k in the Brillouin zone. 

(ii) If at a special point k the Bloch function 9 (r) 
was basis function for a degenerate single-valued 
representation before the perturbation H s was ac¬ 
tivated and if this degeneracy is removed by T-L s , 
then Eq. (7.11) is unusable for the sole reason 
that we do not know which of the basis functions 
of the degenerate representation we should avail 
in this equation. In fact, in this case the Bloch 
spinors i/’k,q,s(i',t) are well defined linear combina¬ 
tions of the functions u s (t)ipk, q {r) comprising all 
the basis functions (/?*,, 9 (r) of the degenerate single¬ 
valued representation (as given, e.g., in Table 6.12 
of Ref. [29]). These specific linear combinations are 
not considered because, at this stage, they are of 
no importance within the nonadiabatic Heisenberg 
model. 

In the framework of the approximation defined by 
Eq. (7.11) the two functions <fik,q,m(r,t) inEq. (7.3) (with 


m = ± 5 ) are usual Bloch functions with the spins lying 
in ±z direction if 

fms(q, k) = S ms - (7.13) 

Otherwise, if the coefficients f ms {q,k) cannot be chosen 
independent of fc, the spin-dependent Wannier functions 
cannot be written as a product of a local function with 
the spin function u s (t) even if the approximation defined 
by Eq. (7.11) is valid. Consequently, even in the limit of 
vanishing spin-orbit coupling, the spin-dependent Wan¬ 
nier functions are neither orthonormal in the local space 
£ nor in the spin space S, but in £ x S only, see Eq. (7.9). 
Thus, also in the case 

Us ->• 0 

spin-dependent Wannier functions clearly differ from 
usual Wannier functions characterized by 

U s = 0 . 

The ansatz (7.5) presents the most general definition 
of Wannier functions. While their localization can be 
understood only in terms of the exact solutions of the 
Schrodinger equation (7.2), the limiting case of vanish¬ 
ing spin-orbit coupling characterized by Eq. (7.11) yields 
fundamental properties of these Wannier functions lead¬ 
ing finally to an understanding of the material properties 
of superconductors [16, 32, 33]. 


VII.2. Symmetry-adapted spin-dependent Wannier 
functions 

We demand that symmetry-adapted spin-dependent 
Wannier functions satisfy in analogy to Eq. (2.15) the 
equation 

1 

2 

P(a)w T ,m(r,t) = Dji(a) ^ d m >m(a)wT',m'{r,t) 

m'=—\ 

(7.14) 

for a £ H since still the assumptions (i) - (iii) of Sec. II. 2 
are valid. Merely the third assumption (iii) is modi¬ 
fied: now the two Wannier functions w T + i(r,t) and 
w T _i(r,f) are situated at the same atom and, conse¬ 
quently, we now put 

WT,m(r,t) = Wi, m {r - R- pi,t), (7.15) 

where m = ±i. 

The vectors T and T' are still given by Eqs. (2.13) 
and (2.16), respectively. The matrices D(a) = [Dij(a)] 
in Eq. (7.14) are again unitary generalized permutation 
matrices, and the subscripts i and j denote the number 
of the atoms at position T and T', respectively, see Def¬ 
inition 2 . 6 . 
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The operators P(a) now act additionally on the spin 
coordinate t, of a function f(r,t), 

P( a )f(r , t) = /(a _1 r - a~H a , a~ l t), (7.16) 

where the effect of a point group operation on the spin 
coordinate t of the spin function u s (t) is given by the 
equation [29] 

u s (a _ 1 f) = d s r s (a)u s r (t) for a £ Hg. (7.17) 

s' 

The matrix 

di/ 2 (a) = [d ss '(«)] (7.18) 

denotes the representative of a in the two-dimensional 
double-valued representation dy/ 2 of 0(3) as listed, e.g., 
in Table 6.1 of Ref. [29], 

We have to take into consideration that the double¬ 
valued representations of a group g are not really rep¬ 
resentations of g but of the abstract “double group” g d 
of order 2\g\, while the single valued representations are 
representations of both g and g d [29]. 

Definition 7.1 (double-valued). Though we use the 
familiar expression “double-valued” representation of a 
group g, we consider the double-valued representations 
as ordinary single-valued representations of the related 
abstract double group g d , denoted by a superscript “d”. 

Since the index in of the spin-dependent Wannier func¬ 
tions is interpreted as spin quantum number, we demand 
that the term 

1 

2 

'y ' d m ' m (a)wT' ,m' (r,t) 

m' = -1 

in Eq. (7.14) describes a rotation or reflection of the 
crystal spin. Thus, we demand that also the matrices 
[dmm' (cc)] are the representatives of the two-dimensional 
double-valued representation dy / 2 , 

Mmm'(«)] = d 1 / 2 (a) for a £ Hg. (7.19) 

Definition 7.2 (symmetry-adapted). We call the spin- 
dependent Wannier functions “symmetry-adapted to the 
double group H d related to space group H ” if they satisfy 
Eq. (7.14) for a £ H d , where the matrices [d mm >(a)] are 
the representatives of the two-dimensional double-valued 
representation dy/ 2 of 0(3). 

Consequently, symmetry-adapted spin-dependent 
Wannier functions are basis functions for the double¬ 
valued representation 

D d = D® dy /2 (7.20) 

of Hg which is the inner Kronecker product of the 
single-valued representation D defined by Eq. (7.14) and 


the double-valued representation dy/ 2 . Thus, the 2p- 
dimensional matrix representatives D ci (a) of D d may be 
written as Kronecker products, 

D d (a) = D(a) x d 1 / 2 (a). (7.21) 

Definition 7.3 (representation defining the spin-depen¬ 
dent Wannier functions). The single-valued representa¬ 
tion D of Hg defined by Eq. (7.14) shall be shortly 
referred to as “the representation defining the spin- 
dependent Wannier functions” and its matrix represen¬ 
tatives 

D(a) = [A»] 

to as “the matrices defining the spin-dependent Wannier 
functions”. 

While usual (spin-independent) Wannier functions are 
basis functions for the representation D defining the 
Wannier functions, spin-dependent Wannier functions 
are basis functions for the double-valued representation 

D d = D ® dy/ 2 

in Eq. (7.20). 

Also the representation D defining the spin-dependent 
Wannier functions has to meet the conditions given in 
Sec. Ill as shall be summarized in 

Theorem 7.1. The two spin-dependent Wannier func¬ 
tion Wy i(r — R — Pi,t) and vu it _i(r — R — pi,t) at the 
position pi are basis functions for the two-dimensional 
representation 

d d = di Cg> d i/2 (7.22) 

of the double group Gg p related to the point group of 
position Gg p . The one-dimensional representations di 
in Eq. (7.22) fix the (generally reducible) representation 
D of Hg defining the spin-dependent Wannier functions 
[Definition 7.3]. The matrix representatives D(a) of D 
still are unitary generalized permutation matrices which 
must be chosen in such a way that they form a represen¬ 
tation of Hg. We again distinguish between the two cases 
(i) and (ii) defined in Theorem 3.2. 

In addition, Theorem 3.3 must be noted. 

Theorem 4.1 does not distinguish between usual and 
spin-dependent Wannier functions but uses only the 
special representations of the Bloch functions or Bloch 
spinors, respectively, at the points k of symmetry. Thus, 
Theorem 4.1 applies to both usual and spin-dependent 
Wannier functions if in the case of spin-dependent Wan¬ 
nier functions we replace the little groups Hy- by the dou¬ 
ble groups H Just as the groups H the groups H^ are 
finite groups in Herrings sense as denoted in Ref. [29] by 
H G' k and, fortunately, are explicitly given in Table 6.13 
ibidem. 

When we consider single-valued representations, then 
the sum on the right-hand side of Eq. (4.1) runs over 
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the p diagonal elements d,(a) of the matrices D^(a) in 
Eq. (5.11). When we consider double-valued representa¬ 
tions, on the other hand, this sum runs over 2p diagonal 
elements df m (a) of the corresponding matrices 

£> d k (a) = S d *(K a )T> d (a)e- zak - t °‘ (7.23) 

where 

S d *(K a ) = S *(K a ) x ( J J ) ( 7 - 24 ) 

[where S (K a ) is given in Eq. (5.8)] because also 
S (K a ) is diagonal and now there are two Wannier func¬ 
tions Wi >m (r — R — Pi,t ) with m = at each position 

Pi- 

We need not to solve Eq. (4.1) directly but we may de¬ 
termine the representations D^ complying with Eq. (4.1) 
in a quicker way. Eq. (7.24) shows that we may write the 

-d 

matrices D fe (a) simply as Kronecker products, 

D fc (a) = D fc (a) x d 1/2 (a), (7.25) 

where D*.(a) is given in Eq. (5.11). 

Now assume that we have already determined accord¬ 
ing to Theorem 4.1 the single-valued representations D]f 
in the closed band under consideration. Then, the rep¬ 
resentations Dk and Djf are equivalent [see Eq. (5.1)] 
and, consequently, also the representations 

D k = Dk®d 1 /2 (7.26) 

and 

D d k = Df®d 1/2 (7.27) 

are equivalent. Hence [Sec. V.l] the double-valued rep¬ 
resentations D ^ comply with Theorem 4.1 in the same 
way as the single-valued representations D]f do. 

Definition 7.4 (affiliated single-valued band). In this 
context we call the band defined by the double-valued rep¬ 
resentations D ^ in Eq. (7.27) the “double-valued band’’ 
and the band defined by the single-valued representations 
D'f an “affiliated single-valued band”. 

While a double-valued band may possess several affili¬ 
ated single-valued bands, any single-valued band is affili¬ 
ated to exactly one double-valued band. 

The affiliated single-valued band is a closed band that, 
generally, does not exist in the band structure of the con¬ 
sidered material. That means that the Bloch functions 
Tk.qi'i") of the closed band under consideration band gen¬ 
erally do not form a basis for the representations 
even if Eq. (7.11) is valid, see, e.g., the single-valued band 
affiliated to the superconducting band [Definition 7.6] of 
niobium as given in Eq. (7.81). 

We may summarize the result of this section in 


Theorem 7.2. Remember that we consider a closed en¬ 
ergy band of p branches and let be given a representation 
D defining the spin-dependent Wannier functions which 
was determined according to Theorem 7.1. The band may 
only be closed after the spin-dependent perturbation H s 
was activated. 

Let be k a point of symmetry in the first domain of 
the Brillouin zone for the considered material and let be 
H^ the little double group of k in Herrings sense. That 
means, is the FINITE group denoted in Ref. [29] by 
H G' k and explicitly given in Table 6.13 ibidem. Further¬ 
more, let be D ^ the 2 p-dimensional representation of 
whose basis functions are the 2p Bloch spinors ifk,q,s{r , t) 
with wave vector k. Deither is irreducible or the di¬ 
rect sum over double-valued irreducible representations of 
H£. The representations D ^ follow Eq. (7.27), 

Dt = Df®d 1/2l (7.28) 

where the p-dimensional representations D)) 1 define the 
affiliated single-valued band. Thus, also each, Dff is the 
direct sum over single-valued irreducible representations 
of H 0 . 

We may choose the coefficients gi q (k ) and fmsid, fe) in 
Eqs. (7-4) and (7.3), respectively, in such a way that the 
spin-dependent Wannier functions are best localized [Def¬ 
inition 2.5] and symmetry-adapted to the double group 
H d [Definition 7.2] if the characters ,\fc( a ) of the single¬ 
valued representations satisfy Eq. (4-1). 

The complex numbers dfia) in Eq. (4-2) stand for 
the elements of the one-dimensional representations di 
of Gq p fixing the given representation D defining the 
spin-dependent Wannier functions [according to Defini¬ 
tion 7.3]. 

VII.3. Time inversion 

VII.3.1. Time-inversion symmetry of the spin-dependent 
Wannier functions 

Within the nonadiabatic Heisenberg model we are not 
interested in spin-dependent Wannier functions that are 
symmetry-adapted to a general magnetic group as given 
in Eq. (6.1), but we only demand that they are adapted 
to the “grey” [29] magnetic group 

M d = H d + I<H d , (7.29) 

or, in brief, we demand that they are adapted to the time- 
inversion symmetry. I\ still denotes the operator of time 
inversion acting on a function of position f(r) according 
to Eq. (6.3) and on Pauli’s spin functions u s (t ) according 
to 

Kufit) = ±u- s (t) (7.30) 

(see, e.g., Table 7.15 of Ref. [29]), where we may define 
the plus to belong to s = and the minus to s = 
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The index in of the spin-dependent Wannier functions 
we still interpret as the quantum number of the crystal 
spin. Consequently, we demand that K acts on m in the 
same way as it act on s, 

Kw T ,m(r , t) = ±w T -m{r , t ) (7.31) 

where again we define the plus to belong to m = +i and 
the minus tom=-|. 

Definition 7.5 (symmetry-adapted to a magnetic 
group). We call the spin-dependent Wannier functions 
“symmetry-adapted to the magnetic group M d ” as given 
in Eg. (7.29) if they are symmetry-adapted to H d [Defi¬ 
nition 7.2], and if, in addition, Eg. (7.31) is satisfied. 

In analogy to Eq. (7.14), Eq. (7.31) may be written as 

1 

2 

Ku’T,m(r,t) = Nii ^2 nm'mWT,m'(r,t ) (7.32) 


and the second Eq. (7.38) is satisfied if n and N in 
Eq. (7.35) follow two conditions, 

d 1 / 2 {a) = nd*/ 2 (a)n _1 f° r a G H d , (7.41) 

and 

D(a) = ND*(a)N _1 for a G H 0 . (7.42) 

The first condition (7.41) is always valid, see, e.g., Ta¬ 
ble 7.15 (q) of Ref. [29], and the second condition (7.42) 
is satisfied if the representation D defining the spin- 
dependent Wannier functions is real. 

In the third Eq. (7.39) the diagonal matrix S ( K) has 
the form 

^\K) = S\K)X (J j) (7.43) 

[cf. Eq.(7.24)] where S *(K) is given in Eq. (5.14). Thus, 
Eq. (7.39) decomposes into two parts, 


where N = [Nij] denotes the ^.-dimensional identity ma¬ 
trix 



(7.44) 


N = 


(1 0 ... 0 \ 

0 1 ... 0 


= 1 


and 

(7.33) 


S *(K) = N*S*(JT)N * _1 (7.45) 


V° 0 ••• 1/ 

and 

n = [rw] = ( J (7.34) 

Eq. (7.32) shows that the 2/^-dimensional matrix 

N d = N x n (7.35) 

is the matrix representative of the operator K of time 
inversion in the co-representation of the magnetic point 
group 

Mq = H([ + KHq (7.36) 

derived from the representation D d in Eq. (7.20). Thus, 
the matrix N d has to comply (Sec. VI) with the three 
equations (6.9), (6.10) and (6.13) which now may be 
written as 

N d N d * = D d (A" 2 ) = - 1 , (7.37) 

D d (a) = N d D d *(a)N d_1 for a G H d , (7.38) 

and 

S d *(K) = N d *S d *(K)N d *-\ (7.39) 

respectively. 

The first Eq. (7.37) is true because 

nn* (= nn) = —1 (7.40) 


which both are evidently satisfied. 

We summarize our results in this Sec. VII.3.1 in 

Theorem 7.3. The coefficients gi q (k) and f ms (q,k) in 
Egs. (7-4) and (7.3), respectively, may he chosen in such 
a way that the spin-dependent Wannier functions are best 
localized [Definition 2.5] and even symmetry-adapted to 
the magnetic group M d in Eg. (7.29) [Definition 7.5] if, 
according to Theorem 7.2, they may be chosen symmetry- 
adapted to H d and if, in addition, the representation D 
defining the spin-dependent Wannier functions used in 
Theorem 7.2 is real. 

Definition 7.6 (superconducting band). If, according 
to Theorem 7.3, the unitary transformation in Eg. (7.5) 
may be chosen in such a way that the spin-dependent 
Wannier functions are best localized and symmetry- 
adapted to the magnetic group M d in Eq. (7.29), we call 
the band under consideration [as defined by the double¬ 
valued representations Df. in Eq. (7.28)] a “supercon¬ 
ducting band v . 

Within the nonadiabatic Heisenberg model, the exis¬ 
tence of a narrow, roughly half-filled superconducting 
band in the band structure of a material is a precondi¬ 
tion for the stability of a superconducting state in this 
material. 

VII. 3.2. Time-inversion symmetry of the matrices f (q,k) 

In this section we derive the time-inversion symmetry 
of the matrices f (q, k) defined in Eq. (7.3) and shall give 
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the result in Theorem 7.4. Thought evidence for this im¬ 
portant theorem was already provided in Ref. [10] and 
later papers [16, 32], we repeat the proof with the nota¬ 
tions used in the present paper. 

Combining Eqs. (7.4) and (7.5) we may write the spin- 
dependent Wannier functions as 

Vh.rn (jP R Pi , t) — 

1 bz n (7.46) 

-4EE e- ik{R+Pi) g iq {k)v k , q ,m{r,i). 

V ^ k g=l 

By application of the operator K of time-inversion on 
Eq. (7.46) we receive 

Kwi >m (r - R pi, t) = 

,bzh (7.47) 

-7= EE eMR+Pi) 9* q {k) K iPk,q, m (r, t). 

fc g=l 

Eq. (7.31), on the other hand, may be written as 
Kw it m(r - R - pi, t) = 

, bz n (7.48) 

T^EE e- ik{R+Pi) g iq (k)v(m)<p k , q ,- m (r,t) 

^ k 5=1 

or, by replacing under the sum k by — k, 


where now the sum runs over —s, and 

+i 

K<Pk, q ,m(r,t) = Y 
s=-§ 

(7.54) 

because [29] 

I<4’k, q ,s(r,t) = v(s)ip- k! q-s(r,t). (7.55) 

Theorem 7.4. Substituting Eqs. (7.53) and (7.54) in 
Eq. (7.52) we obtain the fundamental condition 

fm,(q,-k) = ±f- m ,-e(q,k), (7.56) 

where the plus sign holds for m = s and the minus for 
m = — s. 

Within the nonadiabatic Heisenberg model, the validity 
of this condition is the cause of the formation of sym¬ 
metrized Cooper pairs in superconducting bands [16, 32, 
33]. 

This Eq. (7.56) may evidently be written in the more 
compact form 

i*(q, — k) = nf(g, fc)n _1 (7.57) 

where n is given in Eq. (7.34). 

VII.4. fc-dependence of the matrices f (q,k) 


Kwi, m (r - R- pi,t) = 

BZ ft 

= 7Tf EE e lk{R+Pi) g iq {.-k)v(m)y_k. q ,-m{r, t), 

k 9=1 

(7.49) 

where 

v{±\)=±l. (7.50) 

Comparing Eq. (7.49) with Eq. (7.47) we receive the 

two equations 

9 * q (k) = g iq (-k) (7.51) 

and 

K<Pk, q ,m(r , t) = v(m)<p-k, q ,- m {r, t). (7.52) 

While the first Eq. (7.51) is relatively meaningless, from 
the second Eq. (7.52) we may derive the important 
Eq. (7.56): 

Eq. (7.3) yields the two equations 


Only those bands are of physical relevance in the the¬ 
ory of superconductivity which are closed not before the 
spin-dependent perturbation 7i s is activated. In this sec¬ 
tion we derive the essential property of such bands and 
shall give the result in Theorem 7.5. 

Let be k a point lying on the surface of the first do¬ 
main in the Brillouin zone for the space group H and let 
be H k the little group of k. In this section, k need not 
be a point of symmetry [according to Definition 4.1] but 
also may lie in a line or a plane of symmetry. However, 
we only consider wave vectors k at which Eq. (7.11) is 
valid. Hence, in general, the Bloch functions tp k , q (r) are 
basis functions for a one-dimensional (single-valued) rep¬ 
resentation of Hk- Nevertheless, in very rare cases, the 
Bloch function <p k . q (r) can be a basis function for a de¬ 
generate (single-valued) representation. Both cases shall 
be examined separately. 

Just as in Eq. (3.1) of Ref. [25] we arrange the 2/i Bloch 
spinors u s (t)<p k , q (r) in Eq. (7.11) as column vector 

/ u + i(t)<p k ^(r) \ 
u_i(t)<p k!fi (r) 


+ 5 

<P-k,q,-m(r,t) = Y f-m-s{q,-k)lf- k , q - s (r,t) 
s =-§ 

(7.53) 


$k{r,t) 


u + i(t)ip k , 2 {r) 

u_i(t)<p k , 2 (r) 

u + i(t)ip k ,i(r) 

v u_i(t)ip ktl (r) 


(7.58) 
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with increasing energy, 

Ek,q-1 < Ef~ q < E k , q + 1 . (7.59) 

Then the analogous column vector <f>fc(r,t) consisting 
of the Bloch spinors fpk,i,m ( r , i) in Eq. (7.4) may be writ¬ 
ten as 


$fc( r > t) = g d (k) • f d (fe) • $ fc (r, f) (7.60) 


where 


g d (fc) = g(fe) x(j°) 


(7.61) 


and 


f d (fc) 


/ 0 0 0 \ 

0 0 f (2, k) 0 

\ 0 0 0 f (1, fe) 


(7.62) 


The matrices g (k) and f (q,k) are defined by Eqs. (7.6) 
and (7.7) and still follow Eqs. (2.8) and (7.8), respec¬ 
tively, and 



The matrices g d (k) ■ f d (k) must satisfy Eqs. (4.8) 
and (4.29) of Ref. [25] in order that the Wannier func¬ 
tions are symmetry-adapted and best localized. [We shall 
consider only Eq. (4.29) of Ref. [25] because this equation 
comprises Eq. (4.8) ibidem}. 

Using the notations of the present paper, Eq. (4.29) of 
Ref. [25] may be written as 

Dfc(o) = (g d *(k) ■ f^(fc))- 1 ■ D 1(a) ■ (g d *(k) ■ f d *(k)) 
for a £ (7.64) 


~^~d 

where the matrices D k (a) and D k (a) denote the repre- 

—d 

sentatives of the the representations D k and D k given 
in Eqs. (7.28) and (7.26), respectively. Assume that the 
representations D^, are determined according to Theo¬ 
rem 7.2. Then the representations £)f and D k as well 

—d 

as the representations Dfc and D k are equivalent for the 
points k of symmetry. Consequently, these representa¬ 
tions are even equivalent in any point k of the Brillouin 
zone because the compatibility relations are valid in a 
closed band [25]. First, from the equivalence of £>f and 

D k it follows that the equation 

Df (a) = g*~\k) ■ D fc (o) • g *(k) for a £ H k (7.65) 
is solvable for any k. 


VII.4-1. The Bloch functions ipk,q(r) are basis functions 

for a non-degenerate representation 


In this subsection we assume that the Bloch states 
Vk.qir) are basis functions for a one-dimensional (single¬ 
valued) representation of H k . 

The representations D^. are the direct sum over the 
double-valued representations of the Bloch spinors in the 
considered band, as arranged in the column vector given 
in Eq. (7.58). Hence, the matrices Df a) on the left hand 
side of Eq. (7.64) may be written as 


/ d k fl (a) 0 0 0 \ 


D£(a) 


l 0 


0 

0 


dfc ,2 (a) 0 

0 d fc ,i(a) J 


x d 1/2 (a) 
(7.66) 


(for a = {a\t a } £ H^), where the Bloch state ipk/jir) is 
basis function for the single-valued representations d k q . 

The matrices on the right hand side of Eq. (7.64) may 
be written as 


(g d *(k) ■ f d *(k)) 1 • D 1(a) • (g d *(k) ■ i d *(k)) = 


-rr-d 


(f d *(k)) 1 - [(g rf *(fc)) l -& k (a)-g d *(k)\-t d *(k). 

(7.67) 

Using Eqs. (7.61), (7.25) and (7.65) we may write the 
matrices between the square brackets as 


(g d *(k)) 1 • D fc (o) • g d *(k) 


—d 


(s*~ 1 (k) ■ Dfc(a) • g*(kfj xd 1/2 (a) = 


Df (a) x d 1/2 (a) 

( 0 0 

k.2 ( a ) 


(7.68) 


V 


0 df (a) 0 

aff 


X d 1/2 (a), 


0 d l s x (a) J 


where again the matrices df f/ (a) form single-valued one¬ 
dimensional representations di kq . Remember that [Defi¬ 
nition 7.4] the single-valued representations (T' kq are not 
associated to the Bloch functions of the considered band 
but are fixed by the representation D defining the spin- 
dependent Wannier functions. 

Eq. (7.68) shows that also the matrices between the 
square brackets form a representation being the di¬ 
rect sum over double-valued representations and, hence, 
Eq. (7.64) splits into the n equations 

dk,q ® di/ 2 = k) ■ ® d 1/2 ) ■ f *(q, k ), (7.69) 

(1 < q < n), which are solvable because the represen- 

—d 

tations D J? and D k and, hence, also the representations 
d k ,q ® di /2 and di kq ® d 1 / 2 are equivalent. 
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We now distinguish between two possibilities: 

— If the considered energy band was already closed 
before the spin-dependent perturbation H s was ac¬ 
tivated, then the affiliated single-valued band ac¬ 
tually exists as closed band in the band structure 
of the material under consideration and, thus, the 
representations d kq and d k l q are equal, 

= d% s q . (7.70) 

Hence, all the p, equations (7.69) are solved by 

f(< 7 > fe) = 1, (7-71) 

with the consequence that the Wannier functions 
are, in fact, not spin-dependent but are usual Wan¬ 
nier functions as defined in Eq. (2.6). 

— If the considered energy band was not closed before 
the spin-dependent perturbation 7i s was activated, 
then not all the representations d k q are equal to 
df q . Evidently, the ( 7 th equation is not solved by 
f (q,k) = 1 when d kq ^ d k l q and, consequently, 
the Wannier function actually are spin-dependent. 

We summarize this result in Theorem 7.5. 

Theorem 7.5. If the considered energy band was not 
closed before the spin-dependent perturbation H s was ac¬ 
tivated, the matrices f (q, k) in Eq. (7.3) cannot be chosen 
independent of k. 

In the Sec. VII.5 the matrix f (q, k) shall by determined 
for some points in the Brillouin zone of niobium. 


where d kM and d k l M now are two-dimensional (single¬ 
valued) representations and the matrix f d (fejvf) now is 
four- dimensional, 


f^M) = ( f(l,° :„) ) ' < 7 ' 73 > 

see Eq. (7.62). 

Though d kM ® d 1 /2 and <g> d 1 / 2 again are equiva¬ 
lent, it is not immediately evident that Eq. (7.72) is solv¬ 
able because f d (fcjv/) is not a general 4x4 matrix. How¬ 
ever, also the representations d kM ® d 1 / 2 and d^ M ® d 1 / 2 
have a very special form since they may be written sim¬ 
ply as Kronecker products. Eq. (7.72) indeed is solvable 
since it expresses the most general unitary transforma¬ 
tion between these special representations. 

For instance, consider the point M of one of the bands 
in Table 3 (b) of Ref. [12] and let be d kM = M 3 given 
by the calculated band structure of the material under 
consideration. In addition, let us choose Band 1 in Table 
3 (a) of Ref. [12] as affiliated single-valued band. Thus, 
we have d ^ = M 2 and Eq. (7.72) is solved by 


f d {k M ) 




\ 

0 M 
* 0 ) J 


(7.74) 


as it may be determined by means of the tables given in 
Ref. [29], 

Though both Band 1 and Band 2 in Table 3 (b) 
of Ref. [12] are mathematically correct superconducting 
bands, they cannot be occupied in undoped LaFeAsO [12] 
which, consequently, is not superconducting. 


VII.4-2. The Bloch functions (fik,q{r) are basis functions 
for a degenerate representation 

In rare cases, it can happen that at a special point 
k some of the Bloch states (p k . q {r) are basis functions 
for a degenerate (single-valued) representation and that 
this degeneracy is not removed by the perturbation TL S . 
For example, each of the two superconducting bands in 
the space group P4./nmm = T q D\ h (129) listed in Table 
3 (b) of Ref. [12] consist of two branches degenerate at 
points M and A. The single-valued Bands 1 and 2 in Ta¬ 
ble 3 (a) of Ref. [12] are affiliated to the superconducting 
Band 1 in Table 3 (b) ibidem ; Bands 3 and 4 in Table 3 
(a) are affiliated to Band 2 in Table 3 (b). 

It is crucial for the localization of the spin-dependent 
Wannier functions that also in this case Eq. (7.64) is 
solvable. We reveal the solubility of this equation on the 
example of the bands listed in Table 3 of Ref. [12], 

At point M in each of these bands, Eq. (7.69) may be 
written as 

d kM 0 di /2 = (f d *(k M )y l • «, 0 d 1/2 ) ■ f d *(k M ), 

(7.72) 


VI1.4- 3. Additions 

In this subsection we show that neither Eq. (7.51) nor 
Eq. (7.57) is inconsistent with Eq. (7.64). Remember 
that in this section we only consider points k at which 
Eq. (7.11) is valid. 

First, taking the complex conjugate of Eq. (7.65), we 
receive with D^f* = D^ k and D k = D~ k the condition 

Df k =g-\k)-D_ k - g (k) (7.75) 

showing that we may chose 

g*(-*0=g(fc) (7-76) 

and, hence, Eq. (7.51) is consistent with Eq. (7.75) and, 
consequently, with Eq.(7.64). 

Secondly, transforming the complex conjugate of 
Eq. (7.69) with the matrix n in Eq. (7.34) and using 
d k , q = = d^ q , and 


d\/ 2 — n 1 di/ 2 n 


(7.77) 
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FIG. 1. Band structure of Nb after Mattheis [34]. The dotted 
line denotes the superconducting band. 


(see Eq. (7.41)), we obtain the equation 


respectively. We may unitarily transform the Bloch 
spinors ipk,q,s{r,t) of this single energy band charac¬ 
terized by the representations (7.80) into best localized 
and symmetry-adapted spin-dependent Wannier func¬ 
tions because Theorem 7.2 yields with H 0 = Oh, l-i = 1, 
Pi = 0, Gq p = H 0 = Oh, and d\ = first the single¬ 
valued representations 

Df = T+ Df = H+, Df = P 2 , and D'f = N+ 

(7.81) 

and then, with Eq. (7.28), the double-valued representa¬ 
tions (7.80). 

The representations in Eq. (7.81) define (the only) 
single-valued band affiliated to the superconducting band 
defined by the representations in Eq. (7.80) [Defini¬ 
tion 7.4]. The representation D defining the spin- 
dependent Wannier functions [Definition 7.3] is equal to 

F+ 

1 2 ) 


d—k,q 0 ^ 1/2 — 

(nf(g, Aj)n — 1 ) 1 • (<P s kq 0 d 1/2 ) • (nf(g, fc)n~ 1 ) 

(7.78) 

showing that, in fact, we may chose 

i*{q,-k) = ni(q, fe)n _1 . 

Hence, Eq. (7.57) is consistent with Eq. (7.64). 

VII.5. Example: Band structure of Niobium 

Consider the superconducting band [Definition 7.6] of 
niobium in Fig. 1, as denoted by the dotted line. At the 
four points of symmetry T, H, N , and P in the Brillouin 
zone for the space group 0\ of niobium, this band is 
characterized by the representations 

^25? -^ 25 : and Pi 

of Ol in the familiar notation of Bouckaert, Srnolu- 
chowski and Wigner [26], which may be written as 

r+, H+, N +, and P 5 , (7.79) 

respectively, in the notation of Bradley and Cracknell 
[29] (see Tables 5.7 and 5.8 ibidem) which is consistently 
used in our papers. When we take into account that the 
electrons possess a spin, we receive 


r 5 

0 

di/2 

II 

'-j 

© 

r 8 ' 

Ht 

0 

d 1/2 

= H+ 

© 

Ht 

P 5 

0 

d\/2 

= Pi 

0 

Ps, 

Nt 

0 

di/ 2 

= K- 




Hence, at the points T, H , P, and N the Bloch spinors 
can be transformed in such a way that at each of the four 
points T, H, N , and P two spinors form basis functions 
for the double-valued representations 

(7.80) 


D = T+ (7.82) 

D is one-dimensional since we have one Nb atom in the 
unit cell. The spin-dependent Wannier functions may 
be chosen symmetry-adapted to the magnetic group in 
Eq. (7.29) because Tj is real. 

The Bloch functions of the superconducting band can¬ 
not be unitarily transformed into usual Wannier func¬ 
tions which are best localized and symmetry-adapted to 
0\ since it was not closed before the spin-dependent per¬ 
turbation 1~L S was activated. Thus [Theorem 7.5], we can¬ 
not choose the matrix f(l,fe) in Eq. (7.3) (with q = 1 
since we only have one branch in the superconducting 
band of Nb) independent of k when we demand that 
the Wannier functions are best localized and symmetry- 
adapted. This important statement shall be demon¬ 
strated by an example: 

Consider the point N with the wave vector kjy in the 
first domain of the Brillouin zone for O®. The repre¬ 
sentations d^ l N | and dfe N ,i in Eq. (7.69) are given by 
Eqs. (7.81) and (7.79), 

df N , i = N+ (7.83) 

and 

d hN , i = N+. (7.84) 

Thus, Eq. (7.69) may be written as 

Nf ® di /2 = (f*(l, fcv)) -1 • ® d 1/2 ) ■ f*(l, k N ). 

(7.85) 

This eciuation is solvable since both representations 0 
di /2 and A^ l “ 0 d 1 /2 are equivalent, but it is evidently not 
solved by f (1, fcjv) = 1. In fact, we receive 

f(l,M= J) (7-86) 

by means of Tables 5.7 and 6.1 of Ref. [29]. This is the 
value of f ( 1 , k) also on the planes of symmetry intersect¬ 
ing at N in the neighborhood of N. Further away from 
N, however, f(1, k) may change since it is k dependent. 


T+, H+, P 7 , and N+, 






20 


In the same way, we find 

-\ +i ) (7.87) 

for the points kp on the line F. 

Eqs. (7.86) and (7.87) demonstrate that f (1, k) cannot 
be chosen independent of k in the superconducting band 
of niobium. 


VIII. CONCLUSION 

In the present paper we gave the group theory of best 
localized and symmetry-adapted Wannier functions with 
the expectation that it will be helpful to determine the 
symmetry of the Wannier functions in the band structure 
of any given material. The paper is written in such a way 


that it should be possible to create a computer program 
automating the determination of Wannier functions. 

In this paper we restricted ourselves to Wannier func¬ 
tions that define magnetic or superconducting bands. 
That means that we only considered Wannier functions 
centered at the atomic positions. When other physi¬ 
cal phenomena shall be explored, as, e.g., the metallic 
bound, other Wannier functions may be needed which 
are centered at other positions, e.g., between the atoms. 
It should be noted that Refs. [25], [27] and [30] define 
best localized and symmetry-adapted Wannier functions 
in general terms which may be centered at a variety of po¬ 
sitions pi being different from the positions of the atoms. 
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